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Low temperature dynamics of Ising ferromagnets under finite magnetic fields are studied in terms 
of quantum spin representations of stochastic evolution operators. These are constructed for the 
Glauber dynamic as well as for a modification of this latter, introduced by K. Park et al. in Phys. 
Rev. Lett. 92, 015701 (2004). In both cases the relaxation time after a field quench is evaluated both 
numerically and analytically using the spectrum gap of the corresponding operators. The numerical 
work employs standard recursive techniques following a symmetrization of the evolution operator 
accomplished by a non-unitary spin rotation. The analytical approach uses low temperature limits 
to identify dominant terms in the eigenvalue problem. It is argued that the relaxation times already 
provide a measure of actual nucleation lifetimes under finite fields. The approach is applied to 
square, triangular and honeycomb lattices. 
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I. INTRODUCTION 

Nucleation phenomena are of basic importance in a wide range of metastable systems which typically involve the 
crossing of a free energy barrier that is large compared to thermal fluctuations A Classical examples of such situation 
are the formation of droplets in an undercooled gas or of crystals in an undercooled liquid, whereas numerous analogies 
can also be found in contexts as diverse as, for instance, material science^, astrophysics- and quantum liquids.^ Owing 
to the initial state of these systems, generally produced by a rapid quench from a stable phase, the decay time before 
escaping from metastability may result extremely large at low temperatures. A significant part of the theoretical 
understanding of these relaxation processes has been amply developed in the study of kinetic Ising ferromagnets as 
microscopic (lattice) models of nucleation. In this framework, the metastable phase can be prepared after equilibrating 
the system under an external magnetic field h which is then suddenly reversed. The system therefore evolves towards 
the full minimization of its free energy via the formation of droplets or small clusters of spins aligned with the new 
field direction. These droplets start growing with very small rates until at least one of them exceeds a critical size, i.e. 
a saddle point configuration or a local maximum in the free energy landscape, thus triggering a rapid magnetization 
change in the whole system. This stems in part from the competition between the energy gained by aligning spins 
with the field and the interface energy created in reorienting previously parallel spins, thus escape from metastability 
essentially occurs when the cost of the latter is outweighted by the gain of the former. 

Several analytical studies have been addressed to elucidate the dynamical aspects of these processes in the low 
temperature limils^iSiLSiS while more recently, the actual evaluation of average nucleation lifetimes has been studied 
combining a range of numerical and analytical efforts ^iLiSiiiiiiSiiSi As a further step in this direction, in this work 
we discuss an alternative low temperature procedure (both numerical and analytical), to estimate the relaxation time 
r of Ising ferromagnets evolving through detailed balance stochastic rules. 17 Specifically, we consider both the usual 
Glauber dynamioA& along with a seemingly minor modification of the latter which however yields entirely different 
characteristics at large magnetic fields^i In either case, we construct a quantum spin representation of the evolution 
operator whose spectrum gap (r _1 ) provides a measure of nucleation rates. In line with the general grounds referred 
to above, the implicit assumption allowing for this identification is that the first passage timeii to create randomly a 
critical nucleus is much longer than the characteristic timescale involved in subsequent growth^. Thus, the relaxation 
of the entire system can be expected to coincide with the inverse of the probability of escaping from the metastable 
well. A posteriori, our results will lend further support to this view. 

Another assumption that is usually made in homogeneous systems -and which is crucial for the feasibility of our 
numerical approach- is that multi spin-flip events as well as fusion between subcritical clusters, are vanishingly rare in 
the low-temperature limilji^. This is supported by our analytic work. Therefore, the relevant length scale over which 
the slow part of the dynamic takes place is of the order of a critical droplet size, the first one to nucleate. Although 
on one hand this prevent us from dealing with small field regimes, where the nucleus becomes macroscopic in the 
limit h — ► , on the other hand this enables us to study other /i-regions using numerically accessible clusters, so long 
as the nucleus can be contained in them. This does not presuppose either a precise knowledge of the nucleus size and 
shape (sometimes a conceptual problem of its ownS), or of the most probable path during a nucleation events^, 
so in this regard our numerical and analytical procedures provide a complementary approach to that of absorbing 
Markov chainsii and other related techniques discussed in Refs. [TollTlTIT4lll6j . 
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For two-dimensional lattices and low temperatures T , the average nucleation time t(/i, T) we aim to evaluate has 
been rigorously shown to be parameterizable as 5 - 

T = A(h)exp[pF(h)}, (1) 

thus, the temperature dependence enters solely in the factor = fcgT (hereafter the Boltzmann constant ks is 
set equal to one). Independently of the stochastic dynamic considered, the exponential argument T has been often 
associated with the energy barrier separating the saddle point from the metastable phase. In fact, for the Glauber 
dynamics the analysis of Ref. Q corroborates this issue. However, the results of Refs. [ilHTll ] clearly indicate that 
such interpretation of T is not always correct, even though the geometry of the critical droplets remains unaltered 
by the change of dynamic. In particular, for the modified Glauber (MG) dynamics considered thereSi, and to be 
discussed in the following sections, despite detailed balance the value of T does not actually yield that energy barrier. 
Moreover, under strong magnetic fields the nucleation process led by the MG dynamics turns out to remain active 
(i.e. r > 0), whereas above a maximum field the standard Glauber dynamic just exhibits a fast decay (T = 0). As 
for the amplitudes A of Eq. as well as for those involved in the nucleation times of a variety of system o 1 ' 1 ? , they 
have been usually difficult to evaluate numerically given their subdominant contribution. However, the lack of a fair 
estimation of these prefactors may introduce far reaching theoretical implications in the nucleation picture^. Recent 
efforts have been addressed to remedy this situation in the context of Ising ferromagnets under finite fieldsi2*ii*iii&. 
On departing from the low temperature regime assumed in Q , these amplitudes have revealed a structure of narrow 
peaksi^ which however rapidly collapses towards a piecewise constant function of h in the limit T — > . Their actual 
values also turn out to be dynamic dependenljifiiAi. Here, we numerically estimate the V and A parameters under 
finite fields using both the Glauber and MG dynamics introduced in Ref. 11]. To check the reliability of our numerical 
operational approach we compare its results with those obtained in square lattice s 5 ! 11 , and then proceed further in 
honeycomb and triangular lattices where no results are previously available. One can expect that beside the evolution 
details, the relaxation parameters will be also affected by the lattice structure as it determines the geometry of the 
critical droplets ultimately controlling the nucleation time 5 ^. Apart from the square lattice, their size and shape are 
not known a priori but neither needed in our procedure. Our low temperature analytic work provides results for T 
and A parameters generally in good agreement with their numerical estimates. 

The layout of this work is organized as follows. In Sec. II we recast the master equation governing the probability 
distribution of these processes in terms of a quantum spin analogy whose 'Hamiltonian' provide the appropriate 
transition rates between the original Ising spin configurations. By means of an ulterior non-unitary spin rotation, 
this results in a symmetric representation of the evolution operator. This simplifies considerably the subsequent 
numerical analysis of Sec. Ill in which the spectrum gap of this latter representation is obtained via standard 
recursive techniques^ in several situations. In Sec. IV we develop the analytical approach, in which low temperature 
limits are used to pick out dominant terms in the hierarchy of equations obtained by applying the quantum spin 
Hamiltonian to an appropriate metastable state. We end the paper with Sec. V which contains our conclusions along 
with some remarks on extensions of this work. 



II. DYNAMICS AND OPERATORS 



Let us then consider an Ising ferromagnet with uniform nearest neighbor (NN) interactions J > between the 
spins s = ±1 of a regular (i-dimensional lattice. Under an applied magnetic field h, taken positive from now on, the 
corresponding Hamiltonian reads 

Hi = — J s r s r < — h s r , (2) 

(r,r'> r 

where the first and second sum run respectively over all NN or bond pairs (r, r'), and all spin locations r of the lattice. 
Since H.i actually defines a classical energy functional, the constituents spins do not have a natural dynamics, i.e. 
dt s r = [Ti-i, s r ] = . Consequently, a specific stochastic evolution must be prescribed so as to emulate the interactions 
between the spins and a heat bath, here modeling the fast degrees of freedom not included in the classical Hamiltonian 
Til. As usual, the underlying non-equilibrium dynamics is then approximated by a discrete Markovian process and 
therefore described by a master equation. The latter governs entirely the time evolution of the probabilities P(s,t) 
of finding the system in a certain spin configuration |s) at time t. If W(s — > s') denotes the (time independent) rate 
or transition probability per unit time at which configuration \s) evolves to \s') , the master equation just adopts the 
continuity form 



d t P(s, t) = [ W(s' -> s) P(s', t) - W(s -> s') P(s, t) 

S f 



(3) 
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Because in the context of Sec. I metastability is imposed by an external field, among the decay processes repre- 
sentable by © we restrict our attention to those in which the total magnetization is not preserved^. One of the 
most studied and physically grounded^ examples of this type is the Glauber dynamic^. Its transition rates involve 
Ising configurations differing at most in the state of a spin at a given site r. With the aid of the local field variables, 
which henceforth we define as ip T — ^ Ylir r') Sr ' > these rates can be written as 

W G (s r -> -s r ) = [1 + e 2 ^+ H ^ s '} - 1 , (4) 

where H — h/T . The iteration of these rules eventually bring the system to the Gibbs distribution as they clearly 
satisfy detailed balance in Eq. that is W(s — > s') e~ Hl ^- s > = W(s' — ► s) e~ Wj "i s } . However, other single spin flip 
or Glauber type processes can also be made consistent with these latter condition, so the approach to equilibrium in 
these problems is not unique. As it was referred to above, a recent case of this situation was introduced in Ref. fi"H 
with the aim of clarifying earlier issues of metastable lifetimes. In the MG dynamics proposed there, the effects of 
the J-interactions and the field h are factorized in the transition rates. More specifically, these are given byS£ 

W MG (sr ^ s r ) = [l + e 2 ^ s 'Y 1 [l + e 2 ^-]" 1 . (5) 

Although it can be easily checked that such rates also comply with detailed balance, it will turn out that each of the 
above dynamics behaves quite differently under strong field regimes. 

a. Mean field excursus - Before constructing a more convenient representation for these processes, we pause briefly 
to consider this latter dynamic at a simple mean field level of description. Despite of being quantitatively uncontrolled, 
on the other hand it is able to account for some relevant qualitative features actually occurring in the MG dynamics 
(see Sees. Ill and IV). Thus, in order to decouple the rather involved hierarchy of equations implicit in Eq. we 
simply approximate the local field variables ip r by their mean value z(s)J/T , in turn assumed to be homogeneous. 
Here, ( s ) denotes the average magnetization whereas z stands for the number of NN spins, i.e. the lattice coordination 
number. After inserting the so approximated Wmg rates in Eq. J3J), we readily obtain the magnetization dynamic in 
terms of a non linear differential equation which at finite fields and low temperature regimes reduces to 

d ^)= l+ \Z}° { l J/T i h/T»l. (6) 

Hence, for the region of our interest the relaxation dynamics comes out to be field independent in this scheme. 
Although this is not the actual case below a minimum h- value (see results of Sees. Ill and IV), yet the analysis of 
Eq. 10 pinpoints a genuine difference with respect to the Glauber dynamics. Notice that for this latter, if h/J> z, 
in the limit T — > the master equation results totally decoupled by the Glauber rates (@J, just as if the spins were 
independent. Then, it follows that dt(s) = 1 — (s) , and therefore the time scale of the Glauber problem is of the 
order an elementary step, namely r = I (see also Sec. III). By contrast, in Eq. (6) the magnetization evolves initially 
with a much slower pace as its change is exponentially plunged by the initial metastable phase. In fact, the integration 
of the reciprocal of Eq. © between ( s ) = — 1 and a subsequent magnetization ( s ) > —I involves large escape times. 
More precisely, with the aid of the exponential- integral function Ei ((s)) and its asymptotic expansions^, we obtain 

T ~£j e2zJ/T > t / j<<i - ( ? ) 

This is consistent with a value of T — 2zJ in Eq. Q which later on will be corroborated both numerically (Sees. Ill) 
as well as analytically (Sec. IV). The corresponding amplitudes however result significantly underestimated by this 
mean field simplification which nevertheless, is already able to capture the metastability of MG dynamics, at least 
under strong field conditions. 

b. Quantum spin representations - We now build up an alternative representation of Eq. Q lending itself more 
readily for a numerical study in finite spin clusters which, as pointed out in Sec. I, can embody the nucleation time 
of much larger systems. Firstly, it is useful to recall the matrix elements of the evolution operator Ti associated to a 
generic Markovian process. In terms of transition rates, these elements are constructed asi£ 

(s'\H\s) = - W(s^s') , s^s', (8) 
{s\H\a) = J2 W(s^s'). (9) 

s' 

This permits to think the master equation in imaginary time as a Schrodinger-like representation | P(t) ) = 
e ~ n 1 | -P(O) ) in which the probability distribution | P(t) ) = ^ s P(s, t)\s) evolves according to the action of the 
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evolution operator -here playing the role of the 'Hamiltonian'- on the initial state | P(0) ) (in our case, a metastable 
Gibbs distribution opposing the new field direction). The specific form of TL in either of the above dynamics can 
be straightforwardly found in terms of spin-^ Pauli matrices a and interpreting the local field variables involved in 
Eqs. (@J and JSJ as local field operators tp* 

<r,r'> 

which just for convenience are taken diagonal, say in the cr z -representation. To connect two z-configurations of spins 
differing in the state of site r, and therefore to account for the off diagonal elements ijHJ, we simply project the 
corresponding 'rate operator' (set by in terms of the usual spin raising and lowering projectors o~^~,a~ . For 
example, using the Glauber rates the operational counterpart of (JSJ) will read 

£ \s>) (s>\H G \s) (s\ = - £ { at [1 + e - 2 <*+*>] 1 + a 7 [l + e 2 ^] 1 } . (11) 

s,s' . s^s' r 

Since [<p%, a^] = , the above ordering of application is immaterial. On the other hand, conservation of probability 
requires the emergence of the diagonal elements JJjJ. They basically count the number of ways in which a given 
configuration \s) can evolve to different states \s') through a single spin flip. This can be properly tracked down by 
using the number operators n r = cr^~o~~ along with the weighting of each flip with its corresponding rate. For the 
Glauber case the analog of Eq. © then becomes 

J2\s){s\H G \s){s\ = J2{^[ 1 + ^ +H T 1 + [1-ftr] [l + e- 2 ^)]- 1 } 

s r 

= \ T. [l-^tanh(^+tf)], (12) 

r 

which together with Eq. I jlljl completes the form of TLg- Certainly, the above reasoning is extensible to the MG 
dynamic as well. The related evolution operator TImg of this case thus finally turns out to be 

Hmg= - iscchtf]r[a+ e ff (l + e- 2 ^)- 1 + a' e~ H (l + e 2 ^ 1 " 

r 

+ i £ (l-o; tanh<#) (l-o£ tanh.ff), (13) 

r 

which of course reduces 2 ^ to Tic when h = 0. Also, it can be easily verified that either of these operators remain 
invariant under the spin inversion a z — * — a z along with the field reversal h — ► —h, as they should. Given the rather 
involved manner in which all spins are coupled through the local field operators <f z 7 exact analytic treatments of 
the spectrum of Tic or of H.mg under generic field and temperature conditions may seem unlikely, even in d — l. 28 
However, by exploiting low temperature limits analytic procedures can be developed and applied, as shown in Sec. 
IV. In addition, numerical progress can be made in fair system sizes by means of a suitable similarity transformation 
which we now discuss. 

c. Symmetric representations - As is knownii, the detailed balance property of rates Q and JSJ warrants the 
existence of representations in which TLq and TImg are symmetric and thereby diagonalizable. Moreover, a common 
transformation for that purpose can be found for both dynamics. To this end, we rotate the corresponding operators 
around the z spin direction using a site dependent pure imaginary angle 

cf> r = -i{ip r + H), (14) 

where the ip's are the original scalar fields introduced in Eqs.(@J and This rotation is produced by the non- unitary 
similarity transformation U — e~ % s with S = \ ^ r <f> r a* , which in turn results in the direct product 



u = (V) u T , u r = 



e-H^+H) o 

e i(^+H) 



(15) 



While the diagonal terms of Hg and Hmg remain unaltered by U, it is straightforward to show that 

UafU- 1 =e^' +H Kf. (16) 
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From this latter relation, one can immediately verify that the rotated Glauber operator TC G 
be casted in the symmetric form 



E 



\l-<4 tanh {ifl + H) - of sech (ip* + H) 



UTCgU 1 can finally 



(17) 



whereas the rotated version Ti' MG = UH' MG U 1 of the MG dynamics is also symmetric and comes out to be 



- J2 [ ( 1 _ °f taim fr ) ( 1 — a r taml H ) - a* sech ip z r sech H 



(18) 



The formal analogy with the Schrodinger picture referred to above now becomes more transparent, as all solutions of 
the master equation are necessarily obtained as superpositions of eigenstates with real eigenvalues (or 'energies') 
A > of Hcrmitian 'Hamiltonians'. In particular, the ground states |^o) of both H' G and Tt' MG coincide and are closely 
related to the equilibrium Gibbs distribution. This is because (tp \ = ( s I i s the left steady state of the original 



stochastic operators (notice that their columns add up to zero), and therefore (i/>o| = (ip\U 1 = J2 S ( s \ e 



t—IS \ - I" J 

modulo a normalization factor \/~Z involving the partition function of the Ising energies ©• It is thereby a simple 
matter to check that in our symmetric representation the dynamics of any classical quantity A (which is already 
diagonal in the cr z -representation, such as the magnetization, the energy TLi or any microscopic correlator), can be 
written as 



WW = \ E A & e ~ PHl{ ' } + E e ~ A ' $ \ u ' lA I ^ ) < ^ I u | p(o) ; 



(19) 



A>0 



Thus, we see that the relaxation times discussed throughout Sec. I can be read off from the first excited level of the 
evolution operators constructed so far and whose numerical analysis we next turn to consider. 



III. NUMERICAL RESULTS 



The main advantage of the symmetric representations l|17|) and l|18|) is that their lower eigenmodes, which are just 
the ones dominating the above non-equilibrium terms, can be efficiently computed using recursion-type algorithms 
devised for hermitian matrices, e.g. the Lanczos techniqueSS,. The latter is particularly appropriate to study system 
sizes capable to accommodate critical droplets arising from not too small field regimes. Specifically, for a square lattice 
the critical nucleus is an L x (£ — 1) rectangle of overturned spins gathered to a similar spin on one of its long sides 
of length C = [2J/h] , where [] denotes the integer part£ Hence, in line with the general arguments of Sec. I one 
could expect that for h/ J ^ 0.5 the spectrum gap of at least a 5 x 4 spin system will suffice to yield actual values of 
nucleation times in the low temperature limit. 

Thus, starting from a random initial state but chosen orthogonal to the Gibbs- like distribution \ipo) referred to 
above, we carried out the standard Lanczos procedure in such spin clusters using periodic boundary conditions 
(hereafter, assumed throughout this Section). Let us first consider the Glauber operator l|17|) . In Fig. 1 we show the 
results obtained from its first excitation level Ai, i.e. above equilibrium, when varying the field h/J € (1,4) at low 
temperature regimes T / J ~ 0.2 — 0.4 . The nucleation time parametrization conjectured by Eq. Q here identified 
with 1/Ai, is consistent with both the data collapse in the main panel as well as with the linear behavior evidenced 
in the inset. In particular, the slopes of the latter detect three typical amplitude values which in turn are used as 
scaling factors in the main panel thus producing, as expected, the collapse of different curves. After a least square 
linear fitting of our data, the corresponding relaxation parameters within the above temperature and field ranges are 
found to be 



_ f 16.(1) J — 6.(0) h , 1< j<2, 
~ { 8.(1) J -2.(0) h, 2<^<4, 



0.4(3), 1<^<2, 
A G = I 1.9(4), !} = 2, 

1.3(3), 2<^<4. 



(20) 



It should be mentioned that below T/J ~ 0.1 and h/J ~ 1, the spectrum gap is gradually comparable to the 
numerical propagation of our roundoff errors, while the convergence of the Lanczos recursion becomes slow and 
erratic. Nonetheless, above those regimes, where these problems do not show up, our results are already in fair 
agreement with those of Ref. ^l| as well as with the low temperature analysis of Sec. IV. As conjectured earlier, size 
effects are negligible around this field region, at least judging from 3x3,4x4 and preliminary results in 6 x 4 spin 
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arrays, all of which can enclose the critical droplets depicted in Fig. 1. In this regard, notice that the corresponding 
values of Tq are consistent with both the surface tension and magnetic energy of such droplets, in turn recovering the 
interpretation of T as an energy barrier. Also by approaching the decoupling condition h/J = z from below, the low 
lying levels which were non-degenerate so far, closely approach one another, as they should, whereas Tq — > 0. 

Bolstered by these consistency checks, we now turn our procedure to honeycomb and triangular lattices for which 
these nucleation parameters are not previously available (see also Sec. IV). Due to the roundoff limitations mentioned 
above, we restricted the computations respectively to h/J > 0.5 , T/J > 0.1 and h/J> 1.5 , T/J > 0.3 . For the first 
situation, Fig. 2 displays the results so obtained in an 18-spin honeycomb cluster (schematized by its lower inset). 
These are in line with parametrization l(l)). and for which our numerical estimations yield 

14.(1),- 10.(1),, 0.5<i<l, i$W< «**<!, 

Preliminary tests using 24-spin clusters showed no substantial differences with these results. This conforms with the 
fact that within our accessible lower field bounds, the above Tq 's at most can entail a 5-spin nucleus (assuming 
the usual T- interpretation still holds). However, the shape of such nucleus can not be inferred only from its surface 
tension (14J) as z is not large enough. In contrast, the results of the triangular lattice lend themselves more readily 
for this purpose, at least for the field range shown in Fig. 3. After analyzing 4x4 and 5x4 triangular clusters, in 
this case further cusps in T and amplitude discontinuities are detected, namely 



32.(1) J — 10.(1) h , 1.5 < ^ < 2 , 
20.(0) J- 4.(0) h, 2 < ^ < 4, A G 
12.(0) J- 2.(0) h, 4 < k < 6, 



0. 7(0), 1.5<^<2. 

1. (0), $ = 2, 



0.4(5), 2<^<4, (22) 
1.4(2), §=4, 
1.2(2), 4<^<6. 

Next, we consider the modified Glauber operator l|18fl . In all studied situations, its numerical treatment comes out 
to be numerically more demanding, i.e. spectrum gaps are even smaller than before, particularly below h/J — z — 2 . 
So, we limit our computations to h/J > 1 , T/J > 0.2 for square lattices, h/ J > 0.3 , T/J > 0.2 (honeycomb) and, 
h/J>2, T/J > 0.3 (triangular). Despite these restrictions, the results of Fig. 4 clearly support larger values of T(h) 
than those obtained for the Glauber dynamic. Also, the amplitude values turn out to be different as well as their 
regimes of validity. Specifically, for the square lattice (Fig. 4a), we find 

^-{f ( Z J ' mK ^>l <2, J"*.*'' (23, 

l8 ' (0)J> 7" 2 ' { 1.(0), J>2, 

whereas for honeycomb (Fig. 4b), and triangular (Fig. 4c) systems the respective parameters become 



14.(1) J — 8.(0) h , 



03 <A<i f - 1 ^)' 0.3<A<i 

^ .7 > A i a(a\ h _ i 



mg ~ \, ?J 7 " W ' h " ~ J ' A MG ~ <^ 1.6(6) , § = 1 , (24) 



6.(0) J, &>1 



1.(0), ^>1 



and 

Tmg~\7^\ T:T ' A MG ~ { 1.1(6) , § = 4, (25) 



J 20.(1) J-2.(0)/i, 2<^<4, 
\l2.(0)J, ^>4, 




It is worth remarking that, as before, the square lattice parameters are in reasonable agreement with those of Ref. [llj . 
which lend us further confidence in the identification of 1/Ai with the nucleation time of the system. Since in all 
studied cases Tmg [h) > (h), notice that the usual association of T with an energy barrier no longer applies for 
this dynamicii2*ii Also, these results give evidence that the nucleation process persists at large times and fields, i.e. 
Tmg = 2Jz > 0, as opposed to the Glauber picture where r = 1 beyond h/J — z. Other tests using much larger 
fields suggest an identical behavior (see also Sec. IV), always maintaining a non-degenerate level Ai. 

Finally, we point out that errors throughout all A's might be actually larger than those estimated above, given 
their subdominant contribution to r in Eq. specially within the smaller field regions. This is reflected in the low 
sensitivity of the data collapse to amplitude changes slightly away from their error bands (arising only from upper 
insets). 
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IV. LOW TEMPERATURE ANALYSIS 



Here we develop a low temperature analysis, starting from the quantum formulation of Section II, which provides 
analytic results for the relaxation parameters T and A. The method employs the unsymmetrised Hamiltonian 



{at P-) R + (J2 O + (*r -Pr + )R-(Y, <>) 



(26) 



where rl are the neighbours of r, and = j(l±^), Low temperature versions of the appropriate rates i? ± (m) are 
used. Here, and throughout, to is an integer corresponding to the "total spin" of the neighbours, and ± relates to 
flip up or down. The unsymmetrised form of H. is easier to work with because the low temperature forms of the rates 
there are nicely separated. 

With | n ) the amplitude corresponding to a domain of n up spins in the eigenfunction for eigenvalue s, the eigenvalue 
problem involves a hierarchy of equations relating | n) to | n ± 1), each of which is of the following schematic form: 

[c n R+ {■■■) + d n R- {■■■)- s] \n)=c n R+{---)\n + l)+d n R-{---)\n-l). (27) 

Here the coefficients c n and d n depend on geometric factors of the lattice such as its coordination number z. These 
equations are consistent with an equal amplitude eigenstate with eigenvalue s — 0. We want the next eigenvalue, i.e. 
the "gap" s = (Ae^)" 1 . 

The relaxation from a metastable state, which we take to have all spins down (i.e. antiparallel to the field), is 
governed by the slow rates, especially the slow up-flip rates. Which rates are small depends on the field, so different 
field regimes have to be considered separately. The first equation (n = 0) has do = 0, i.e. the only terms come from 
cqR + (■■■) = R + (—z), corresponding to nucleation of a single up spin, which state has amplitude |1). The only case 
where R + {—z) is not small at low temperatures is the Glauber dynamics case with zJ < h, where R + (m) ~ 1, all to. 
The equations then give the gap of order 1, so 

T G = , h/J > z. (28) 

In all other cases R + (—z) is exponentially small in (3 at low temperatures, and this results in non-zero T for 
Glauber dynamics (G) with z J > h and for modified Glauber dynamics (MG) at any h. So we confine our attention 
hereafter to those cases, at very low temperatures. Then, R + (—z) is by far the smallest of the flip-up rates (since 
for any positive integer I, for MG dynamics R + (—1) ~ e - 2 0l J <g; \ t while in the Glauber case, if h < I J then 
R + (-l) ~ e -mu-h) < Y y Consequently we may neglect terms involving further factors of R + (—z), as occur, 
corresponding to further nucleations of isolated single-spin clusters, in the equations for n > 0. As a result, in 
Eq. H27fl . for amplitudes |n > 2), it is only necessary to consider "connected clusters" where all up-spins have at 
least one up-spin neighbour. For example the second equation (n = 1) has CiR + (---) = zR + { — (z — 2)) and 
diR~(- ■ ■) = R~(—z) after neglecting the further nucleation terms involving R + (—z). 

While basic ideas and procedures are similar for G and MG dynamics, because of the different forms of their 
rates, the ordering of terms in the equations can be different in some regimes, so we discuss the two cases separately, 
beginning with the Glauber case. 



A. Glauber dynamics 

Case (a): d-dimensional lattices with zJ > h> (z — 2)J. Here R + (— z) ~ e -W(zJ-h) ^ ^ anc j a n other R + (m) ~ 1 
while R~{-z) - 1. 

The resulting recursion type eigen-equations have (as always) the equal amplitude solution with s — 0, and (because 
the only small R + occurring is in the first equation) the next eigenvalue satisfies (in the low temperature limit) 
(R + (-z) - s)(z + l - s) = R+(-z), giving s = z {z + i)-i e -2/3(*J-h) = (Ae^y 1 . So 

T G =2{zJ-h) , Ag = (z + 1)/z, (29) 

for zJ > h > (z — 2) J in any lattice. (This is consistent with the numerical predictions, and it includes the linear 
chain result T G = AJ-2h,A G = 3/2)M 

Case (b): d > 1-dimensional lattices in the next regime, (z — 2) J > h. 

Now, as well as R + (—z), a second flip-up rate, R + ( — (z — 2)) becomes very small. The geometry of the domains 
of n and n± 1 sites determines the numbers to in the rates R + (m) occurring in the equation for \n). In particular the 
size, q, of the smallest ring of bonds on the lattice determines as q— 2 the number of successive equations in which the 
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only up rate is — (z — 2)). The consequence for the gap is that s is proportional to R + (— z) [R + ( — (z — 2)) } q 2 , 
giving the result: 

r G = 2 (z J - ft) + 2 (q - 2) [ (z - 2) J - ft] , (30) 

where q = 4, 6, 3 for square, honeycomb, and triangular lattices. These results agree with the numerical ones. 

The determination of A's is most easily carried out by considering the (first order, nonlinear) recurrence relations 
for n n — \n)/\n — 1). In the low temperature limit, ratios of small rates make the deciding /z n 's tend to 1, for the 
(q - 2) lowest n's, and to for the next one; and in the limit the ratios of the coefficients c n ,d n in that range of 
n's determine the numerical factor in (1 — ^i) and hence Ac- The (analytic, numerical) results are Aq = (3/8, 
0.4(3)) , (l/6, 0.1(6)) , (l/3, 0.4(5)) for square, honeycomb, and triangular lattices, respectively; the agreement is 
good except for the last one. 

The boundaries of the region are set by where new combinations of rates become limiting. The analytically 
determined ranges of validity of the results in this case (b) are 4J > ft > 2 J for triangular, 2J > ft > J for square, 
and J > ft > J/2 for honeycomb lattice. 

Case (c): For the triangular lattice there is further regime (2 J > ft) where yet another up-flip rate, namely 
R + ( — (z — 4)) , becomes small. 

Here we expect, in analogy to the argument and results given above for the previous regime that s will be proportional 
to R + (— z) [ R + ( — (z — 2)) ] 9 2 [ R + ( — (z — 4)) 1 ^ with Q an integer related to topological features of the triangular 
lattice. This gives the form 

T G = 2(z J- ft) + 2(q- 2) [(z - 2) J -ft] +2Q [(z - 4)J - ft] , (31) 

(with z = 6, q = 3 ). The numerical results are consistent with this with Q = 3. 



B. Modified Glauber dynamics 

Procedures for the MG dynamics are in principle similar. But now the primitive rate R + (—z) ~ e -^P^J j s always 
small, in all regimes (even h very large), so always 



MG 



* 0. (32) 



For the following we use the notation e = R+ (-z); A = R+ ( - (z - 2)); 7 = R+ ( - (z - 4)); 6 = R~ (-z) ~ 
i?~ ( — (z — 2)); yu„ = \ n)/\ n— 1) = (1 — A„). Then for any lattice the first two recurrence equations are 

s = e\i, (33) 
zA\ 2 = SXx/fix + s. (34) 

After these, the equations become lattice-dependent; e.g. for the triangular lattice the next ones are 

2 7 A 3 = 2(5A 2 /M2 + s, (35) 
3 7 A 4 = 3SX 3 /fi 3 + Sl (36) 
47A5 = 3(5A 4 /M4 + s, (37) 

etc. Details of the further reductions depend on the field regime. 
Case (a): For ft > (z — 2) J, the rates satisfy 8 -c A <C 7. 

Then for all the lattices the equations give two possible forms of solution: A„ <gi 1, so fj, n <~ 1 and s — (ground 
state); or fi n <C 1 , so s = e ~ e _2/3z,/ . Hence 

Y MG = 2zJ , A MG = 1. (38) 

This applies for any lattice, including the linear chain. 

Case (b): In the next regime (z — 4) J < ft < (z — 2) J, possible in d > 1-dimensional lattices, the rate ordering is 
A < 5 < 7. 

First consider the specific case of the triangular lattice. Again because S <C 7, for s ^ the higher equations of the 
hierarchy give A„ <~ l,n = 2,3, ... But now A <C 5, so the second equation of the hierarchy gives Ai = zA/S and it 
follows that s — zeA/S. Hence 

Amg = 1/z, r M G = 2Jz + 2J(z-2)-2ft = 20J-2ft, (39) 
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for the triangular lattice in this regime. 

For the other lattices, a cycle of (q — 2) successive equations (after the first equation) involve A/S and that is the 
origin of the general form 

r MG = 2zJ + 2(q-2)[(z-2)J-h]. (40) 

The (analytic, numerical) results for Amg f° r {z — 2)J > h, arc (1/8, 0.2(3)) , (1/6, 0.1(3)) , (1/6, 0.3(4)) for square, 
honeycomb, and triangular lattices, respectively. 

Case (c): On the boundary h = (z — 2) J between the last two regimes. Here the (analytic, numerical) results for 
Amg are (11/8, 1.4(2)) , (11/6, 1.6(6)) , (7/6, 1.1(6)) for square, honeycomb, and triangular lattices, respectively. 

V. CONCLUSIONS 

The low temperature relaxational kinetics of Ising ferromagnets in a field has been treated for various lattices for 
both Glauber(G) and modified Glauber (MG) processes using a quantum representation of the evolution operators. 
The unsymmetrised version is convenient for the analytical work in the low temperature limit (Section IV) while the 
recursive numerical approach (Section III) requires the symmetrised form produced by a special spin rotation. 

The gross features produced by the two approaches agree completely. Those features include striking differences 
between the Glauber and modified Glauber cases: the latter case is always activated, even at high fields; in each of a 
sequence of field regimes, for both processes T is found to be a linear function of h and the amplitude A is a constant, 
but both A and the form of T, and also the field regimes, differ between Glauber and modified Glauber cases. 

While A and the slope of T(h) are discontinuous at regime boundaries, T is continuous. This can be understood 
from the analytic discussion, where it is seen that the regimes are distinguished by which rates are most limiting, 
and in the low temperature limit the exponents in the limiting rates cross over at the regime boundaries, and those 
exponents determine V but not A. 

Both the numerical and analytic investigations make no use of the shape of critical droplets nor of the most probable 
path towards a nucleation event, and were carried out for square, triangular, and honeycomb lattices. The numerical 
approach indicates, and the analytic one confirms, that T is lattice-dependent (except in the non-activated high field 
regime of the Glauber case, where it is zero), and also A is lattice-dependent except in the highest field regime for the 
modified Glauber case. According to the analytic work, for the highest field activated regime the lattice-dependence 
of r involves just the coordination number z for both Glauber and modified Glauber cases, but in subsequent regimes 
other geometrical aspects of the lattice, e.g. the smallest ring size q, affect the value of T. Similarly the lattice- 
dependence of A is, as one moves down the field regimes, first through z, and thereafter involving further aspects of 
the lattice. 

The quantitative agreement between the predictions of the numerical and analytic approaches is very good for T, 
and for the regime ranges, and slightly less good for A, particularly for the smallest field regimes. This is as might 
be expected, since (i) at low temperatures a given (but bounded) numerical error in the gap evaluation will mostly 
propagate an error in A, rather than in T, as the latter carries an extra weight proportional to (3; and (ii) as h 
decreases the gap becomes smaller, and so does the accuracy of the machine calculations. To understand more fully 
the comparisons between the numerical and analytic results it would be desirable to generalise the analytic work to 
finite low temperatures. 

The analytic work includes predictions for arbitrary lattices (e.g. for the regimes in which the results for T, A 
depend only on z). It would be valuable to extend this, and to extend the numerical work to other, especially three- 
dimensional, lattices. Further suggested extensions of the work are to other models (e.g. Potts models, where domain 
walls remain sharp) and to disordered cases: even a low concentration of weak bonds can make nucleation much faster. 
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FIG. 1: Low temperature estimation of the relaxation parameters F = T\xi{t/A) in square lattices for T — 0.2 (squares), 0.3 
(triangles) and 0.4 (circles) resulting from the Lanczos diagonalization of Eq. I17H in small clusters (up to 6 x 4 spins). Solid 
lines denote the estimations of Eq. (121 )t . The data collapse was attained upon using the amplitudes A — e s derived from the 
slopes s of the upper inset. From top to bottom they refer respectively to h = 1.7, 2 and 2.3 , characterizing typical regimes of 
1171 . The lower lines sketches the shape of the critical droplets for 1 < h < 2 (three spins), and 2 < h < 4 (single spin). 
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FIG. 2: Relaxation parameters of honeycomb lattices estimated from the diagonalization of Eq. 1171 using the 18-spin cluster 
depicted by the lower inset. As in Fig. 1, the amplitudes yielding the data collapse in the main panel [T = 0.2 (squares), 0.3 
(triangles) and 0.4 (circles) ], were inferred from the slopes of the upper inset. The latter refer respectively to h = 0.75 (top), 
1 (middle), 1.5 (bottom) and are representative of the field regimes summarized in Eq. Ill' 1 1 (solid lines of main panel). 
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FIG. 3: Relaxation parameters of triangular lattices arising from the gap of Eq. 1171 in 4 x 4 and 5x4 spin clusters. As before, 
the data collapse was obtained from the slopes of the upper inset. The former refers to [T = 0.3 (triangles), 0.4 (squares) and 
0.5 (circles) ], and follows closely the field regimes given in Eq. (1221 . denoted by solid lines. The amplitudes resulting from the 
inset slopes are characteristic of the regimes identified in Eq. 1221 . Here, they refer to h — 1.8, 2, 3, 4, and 4.5, in descending 
order. The size and shape of critical droplets are schematized below. From left to right they refer to 5-spins (1.5 < h < 2), 
2-spins (2 < h < 4) and a single spin (4 < h < 6). 
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FIG. 4: Relaxation parameters of (a) square, (b) honeycomb and, (c) triangular lattices under the modified Glauber dynamics 
given in Eq. @. Solid lines in each case stand respectively for the regimes identified in Eqs. 1231 . 1241 and, 1251 . In (a) and (b) 
they follow closely the data of T = 0.4 (circles), 0.3 (triangles) and, 0.2 (squares), obtained from the numerical diagonalization 
of Eq. 1181 . Similarly for (c), where the points denote T = 0.8 (triangles), 0.6 (circles) and, 0.4 (squares). In contrast to the 
standard Glauber dynamic, here the relaxation process remains active at large fields, i.e. T > for h 3> J ■ As before, the 
amplitudes were derived from the slopes of the insets. In descending order they refer respectively to (a) h = 1.5, 2, 2.5 , (b) 
h = 0.75, 1 , 1.5 and (c) h = 2, 2.5, 4, 5 (horizontal line), typical cases of each situation. 



